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The radiation background produced by the 21 cm spin- flip transition of neutral hydrogen at 
high redshifts can be a pristine probe of fundamental physics and cosmology. M, z ^ 30-300, the 
intergalactic medium (IGM) is visible in 21 cm absorption against the cosmic microwave background 
(CMB), with a strength that depends on the thermal (and ionization) history of the IGM. Here we 
examine the constraints this background can place on dark matter decay and annihilation, which 
could heat and ionize the IGM through the production of high-energy particles. Using a simple 
model for dark matter decay, we show that, if the decay energy is immediately injected into the 
IGM, the 21 cm background can detect energy injection rates > 10~^* eV cm"'' sec~^. If all the 
dark matter is subject to decay, this allows us to constrain dark matter lifetimes < 10^^ sec. Such 
energy injection rates are much smaller than those typically probed by the CMB power spectra. The 
expected brightness temperature fluctuations at z ~ 50 are a fraction of a mK and can vary from the 
standard calculation by up to an order of magnitude, although the difference can be significantly 
smaller if some of the decay products free stream to lower redshifts. For self-annihilating dark 
matter, the fluctuation amplitude can differ by a factor < 2 from the standard calculation at 
z ~ 50. Note also that, in contrast to the CMB, the 21 cm probe is sensitive to both the ionization 
fraction and the IGM temperature, in principle allowing better constraints on the decay process and 
heating history. We also show that strong IGM heating and ionization can lead to an enhanced H2 
abundance, which may affect the earliest generations of stars and galaxies. 

PACS numbers: 95.35.-|-d, 98.62.Ra, 98.70.Vc 



I. INTRODUCTION 

The cosmic "dark ages," stretching from the last 
scattering surface of the cosmic microwave background 
(CMB) dX z ^ 1100 to the formation of the first lumi- 
nous sources at z ~ 30, are one of the last frontiers for 
observational cosmology. The problem is the difficulty 
of finding any probes: because the intergalactic medium 
(IGM) has mostly decoupled from the CMB, there is no 
background radiation field against which we can study 
it, and (by definition) there are no local light sources. 

This is unfortunate because, at least in principle, the 
physics of the dark ages is sufficiently simple that we 
can hope to understand it in detail. The only factors 
that enter are the CMB, the expanding Universe, recom- 
binations, Compton scattering (which couples the CMB 
to the IGM), and gravitational growth - almost entirely 
in the linear regime. Thus, the IGM properties during 
the dark ages (and especially its fluctuations) constitute 
a strong test of fundamental cosmological parameters in 
an entirely analogous way to the CMB. Conversely, if we 
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take the physics as well-understood, probes of the dark 
ages would offer stringent tests of exotic physics. 

Two examples are dark matter decay and annihilation. 
If either of these processes happen (even for only a frac- 
tion of the dark matter), they could inject high-energy 
photons (or other particles) into the IGM. These would 
then scatter and deposit some or all of their energy as ion- 
izations and heat. An altered reionization history may 
affect the total optical depth and, as a consequence, the 
CMB power spectrum. Decaying dark matter particles 
were in fact initially advocated as a possible explanation 
for the high optical depth detected in the first year of 
WMAP data [ElllllllEll- Such data allowed an up- 
per limit on the energy injection rate and consequently 
on the decay time {t > Hq^). (These constraints, and 
others mentioned throughout this paper, were obtained 
with the first-year WMAP data, but we do not expect the 
newer data ,7| to improve them substantially given the 
many uncertainties on the various reionization processes 
affecting the total optical depth.) 

However, the CMB is only affected if a substantial frac- 
tion of the CMB photons interact with the IGM - or in 
other words if the optical depth to electron scattering 
is substantially altered by the decay products 0, H, 0] ■ 
Compared to the IGM temperatures of ~ 10-10^ K at 
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z ^ 30-300, ionization requires a substantial energy in- 
put, and decays with long timescales (or low energy in- 
jection rates) cannot be ruled out with the CMB. Never- 
theless, because the IGM temperature is so small, these 
models can still significantly affect the thermal history of 
the IGM and hence the history of star formation. A 
number of recent particle physics models motivate such 
scenarios (see the discussion below in Ap . 

Dark matter annihilation has also been considered 
recently, primarily as an explanation for gamma-ray 
sources in the local Universe (e.g., 9]). But any such sce- 
nario will also predict an annihilation background from 
encounters between IGM particles; because the physical 
density increases like (1 -I- zf" , such events could be rel- 
atively common in the early Universe and so affect the 
CMB JUJ. Of course, these scenarios will continue to 
affect the IGM throughout the dark ages as well. 

The 21 cm background can offer powerful constraints 
on dark matter decay and annihilation (or indeed any 
exotic process that injects energy into the IGM during 
the dark ages). The hyperfine level populations of the 
ground state of neutral hydrogen are determined through 
competition between absorption of (and emission stimu- 
lated by) CMB photons (and possibly UV photons) and 
collisions. Since the IGM is dense at high-rcdshifts, col- 
lisions dominate, driving the level populations into equi- 
librium with the kinetic temperature of the gas. Because 
the latter is colder than the CMB, the IGM is a net ab- 
sorber of CMB photons at z - 30-300. Thus, the 21 cm 
transition can be used to map ffuctuations in the IGM 
during the dark ages, offering a pristine probe of cosmol- 
ogy 11]. Previous work has emphasized the possibility 
of constrai ning the matter power spectrum with this tool 
[Til IT^ Il4| . Here we point out that the fundamental 
properties of the fluctuations - on all scales - depend sen- 
sitively on the thermal history (see also Ref. 15). Thus 
the 21 cm background can be used to constrain dark mat- 
ter decay and annihilation. It is much more sensitive than 
the CMB to dark matter decay because (i) it depends 
on the thermal history, not just the ionized fraction, and 
(ii) it is directly sensitive to the late time behavior (when 
most of the energy is injected). It is less useful for con- 
straining annihilation scenarios because in that case a 
large fraction of the energy is injected at early times. 

The excess heating and ionization induced by dark 
matter decay and annihilation can also affect the chem- 
istry of the IGM. In the simple chemical environment of 
the primordial IGM, their most important effect will be 
on the II2 abundance, which is significant because that 
molecule is an important coolant for low-temperature gas 
[3)01 and is th ough t to pla y a key role in the formation 
of the first stars [H [H [M Ml ■ Unfortunately, the im- 
plications - and even magnitude - of any possible boost 
to early structure formation are controversial 10, 22, 
Here, we provide a more detailed look at H2 formation 
and incorporate some hitherto neglected effects in the 
calculation. 

The rest of this paper is organized as follows. We 



briefly review the physics of the 21 cm background in 
01 We describe our simple model for dark matter decay 
and annihilation, and review some particle physics mo- 
tivation, in ijlHI We then present our results for the 21 
cm background in illVI and for the H2 abundance in f0 
Finally, we discuss their implications in MVll 

In our numerical calculations, we assume a cosmol- 
ogy with f7„ = 0.26, f^A = 0.74, = 0.044, H = 
100/ikms-i Mpc-i (with h = 0.74), n = 0.95, and 
(Tg = 0.8, consistent with the most recent measure- 
ments Q , although we have increased cts from the best- fit 
WMAP value to improve agreement with weak lensing. 
We quote all distances in comoving units, unless other- 
wise specified. 



II. THE 21 CM BACKGROUND 

We review the relevant characteristics of the 21 cm 
transition here; we refer the interested reader to Ref. [^ 
(and references therein) for a more comprehensive dis- 
cussion. The 21 cm brightness temperature (relative to 
the CMB) of a patch of the IGM is 
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where 5 is the fractional overdensity, x^a = \ — Xi is the 
neutral fraction, Xi is the ionized fraction, Ts is the spin 
temperature, T^ is the CMB temperature, and du||/dr|| is 
the gradient of the proper velocity along the line of sight. 
When Ts < T^, the IGM appears in absorption. The last 
factor accounts for redshift-space distortions |T^l25| . 

Before the first luminous sources turn on, the spin tem- 
perature Ts is determined by competition between scat- 
tering of CMB photons, collisions [2^, and scattering 
of Lya photons [53, [2^ [l^- In equilibrium (which is 
achieved rapidly). 
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Here Xc is the total coUisional coupling coefficient, in- 
cluding both H-H interactions [13, [111 and H-e~ colli- 
sions 13^ I33I I . We will denote the separate coefficients 
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and 
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respectively. The last term describes 



Wouthuysen-Field coupling: Xa is the coupling coefficient 
and Tg is the effective color temperature of the radiation 
field [2i[33]. Typically T^ w Tk in the IGM, and 
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where Sa is a factor of order unity describing the detailed 
scattering process [2^, [3i [IM 36] and J a is the radiation 
background at the Lya frequency, in units of photons 
cm~^ s" 
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From equation it is obvious that fluctuations in the 
density, temperature, ionized fraction, radiation back- 
ground, and velocity all source fluctuations in the bright- 
ness temperature. Because, to linear order in fc-space, 
velocity perturbations are simply proportional to density 
perturbations, we can write the Fourier transform of the 
fractional 21 cm brightness temperature perturbation as 

J2i(k) = ip + ^i^)6 + (3h6h + PaSa + PtSt, (4) 

where /i is the cosine of the angle between the line of 
sight and the wavevector k and each Si describes the 
fractional variation in a particular quantity: da for the 
Lya coupling coefficient Xq,, Sh for the neutral fraction, 
and St for Tk- The expansion coefficients (3i are 
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Tk - Xtot{l + 2:tot) dlnTft- ' 

where xtot =Xc + Xa- 

III. DARK MATTER AND THE IGM 

The implications of dark matter decay (and/or an- 
nihilation) for the thermal and ionization histories of 
the IGM have been considered by a number of authors 
P, 0, 0, 0, H, S Q • Here we first review some particle 
physics scenarios that can affect the dark ages, and then 
we show how to compute the resulting IGM histories. 

A. Some Example Particle Physics Scenarios 



from the standard case hy z ~ 100 and at z ~ 20 is ~ 20 
times higher, if energy transfer is perfectly efficient 0. 
However, in these particular models a detailed account- 
ing of the interactions between decay products and the 
IGM substantially decreases their effect and pushes it to 
lower redshifts j23 |. 

Another possible candidate is a sterile neutrino with 
rrixc^ ^ 2-4 keV. Again, even for the highest mass in 
this range, the ionized fraction produced by these parti- 
cles is small (x^ < 0.1 at z = 0), but the present-day IGM 
temperature increases by several orders of magnitude (al- 
though again the energy transfer is likely not perfectly 
efficient; 0, ^^)- Although such a neutrino is unlikely 
to be the sole dark matter component js^, it may still 
constitute a non-negligible fraction of the dark matter. 
Such neutrinos have decay times tx —3 x 10^^ sec. 

Decays of super-heavy dark matter particles {mxc^ > 
10^^ GeV) have been considered as possible sources of 
reionization |39l |40| and also as a production mecha- 
nism for ultra-high energy cosmic rays With 
such energetic particles, the precise ionization and heat- 
ing strongly depend on the type of reactions assumed. 

Two kinds of annihilating dark matter have been in- 
vestigated in connection with non-standard reionization: 
neutralino annihilation (e.g. [i^li^ l'). involving particles 
with mxc^ > 30 GeV, and light dark matter annihilation 
(mxc^ ^ 1-100 MeV ^]). While the expected ioniza- 
tion fraction is small, it extends to high redshifts, greatly 
affecting the CMB power spectrum. These models can 
already be constrained with current CMB data and will 
be better constrained by future CMB polarization mea- 
surements 0, 01 . In the most reasonable models, the 
IGM temperature can be several times larger than nor- 
mally expected at z ~ 20-30. 

The effects of light dark matter annihilation on CMB 
power spectra have also been considered. CMB data still 
allow models with small ionization fractions but with 
IGM temperatures a few times larger than in the stan- 
dard scenario (e.g., ~ 2 times larger at z = 60; 



Dark matter particles can ionize and heat the IGM 
through either decay or annihilation. 

For decay, short hfetimes {tx < 3 x 10^® sec, with X 
denoting the dark matter particle) are disfavored, as they 
cause the CMB power spectra to disagree with exist- 
ing data P, 1^. On the other hand, particles with ex- 
ceptionally long lifetimes, such as gravitinos with mass 
mxc^ ^ 10-100 MeV (and hence tx ^ lO^^ sec [13) only 
affect the IGM temperature at low redshift (z < 2) 
We wish to consider particles with lifetimes in the range 
tx - 1024-1027 sec. 

One such example is an axino with mass mxc^ ^ 1- 
100 MeV; positrons produced as their decay products 
could explain the 511 keV emission excess from the 
Galactic center |33- Axinos in this mass range have 
ijf ~ 3 X 1024-1026(MeV/mxc2)sec. The ionized frac- 
tion produced by such particles is too low to be observ- 
able with the CMB 8] , but the IGM temperature departs 



B. Energy Deposition 

Rather than examine each of these models individually, 
we will follow the method of Ref. who presented a 
simple but general model for the effects of dark matter 
decay. We first suppose that the dark matter particle 
has a decay rate Tx = ^x^i ^'^ that the physical number 
density of particles is 

nx{z)^n\{l + zfe'^-\ (9) 

where rfi^ is the comoving number density. In our calcu- 
lations we will assume Tx «C Hq, so that the exponential 
factor can be ignored. While not strictly necessary, this 
is probably the most interesting case (because otherwise 
the decaying particle could not serve as the dark matter 
we see around us today) and serves to illustrate the basic 
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results. The total comoving emissivity is then 

Q^^^'^^ = frxnxmxc\ (10) 

where / is the fraction of decay energy that is potentiaUy 
available to the IGM (e.g., excluding particles that free 
stream without interacting, such as neutrinos) ; the mean- 
ing of the superscript will become apparent momentar- 
ily. It is convenient to define £,x = ^x^x/^b ~ 5.9Tx 
to be the decay rate normalized for notational conve- 
nience to the baryon fraction (i.e., the equivalent de- 
cay rate if it were the baryons that decayed). Then 
Qthick _ f£,xn^iTipC^ , where is the comoving number 
density of baryons. 

For concreteness, we will assume that the particle de- 
cays to high-energy photons, which then either scat- 
ter through the IGM or free stream to lower redshifts. 
(Other decay products are of course possible and will 
have similar qualitative effects; the important parameter 
is the total energy injection rate.) High-energy photons 
can interact with their surroundings through photoion- 
ization, Compton scattering, pair production (via scat- 
tering off neutral atoms, free electrons, CMB photons, 
or ions), and scattering with CMB photons the 
high redshifts of interest here, photoionization and scat- 
tering off of free electrons and atoms dominates energy 
loss at low photon energies {E < 3 keV); in this regime, 
the IGM is opaque. At high energies (E > 10-10^ GeV, 
depending on redshift), pair production off CMB photons 
dominates, and the total optical depth is also large. The 
resulting pairs rapidly deposit their energy in the IGM. 
However, between these thresholds lies a "transparency 
window" of relatively low optical depth; photons in this 
regime experience optical depths (over a Hubble length) 
T ^ 10~^-1 at z = 100 1]. Thus they lose only a fraction 
of their energy to the IGM, with the rest free streaming 
to the present or being lost to cosmic expansion. 

The decay energy can be deposited in heating, coUi- 
sional excitation, or ionization; while the fractions Xh, 
Xe , and Xi in each process actually depend on the initial 
photon energy and the ionized fraction |48j |. to the accu- 
racy of our simple models the rules Xi ~ Xe ~ — Xi)/3 
and Xh ^ + 2xe)/5 suffice 0]. 

When the decay products experience a large optical 
depth, the energy deposition rate is simply Q*''"^"^ (hence 
the superscript). However, if the energy is injected in the 
transparency window, we have instead 
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where r(z) = cn},{z)aes/ H{z) is the effective opti- 
cal depth for interactions over a Hubble length. Here 
(Jeff is the energy-averaged cross-section for interactions. 
Rather than attempt to compute the relevant processes 
in detail, we will write 
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where rioo is the optical depth at 1 -I- z = 100 and the 
exponent assumes a constant doff (so that r cx Ub/H). 



More detailed calculations of r(z) can take into account 
the specific interactions with the IGM for a given par- 
ticle physics model (e.g., |23)j but we forego these here 
to keep our calculations generic. Thus decay into the 
transparency window can be described in the same way 
as instantaneous energy injection provided we make the 
replacement S^x — > S,x[^ — exp(— r)]. In this case, most 
of the energy would free-stream to the present day and 
contribute to the X-ray and 7-ray backgrounds. If the de- 
cay photons are monoenergetic, the observed background 
provides fairly powerful constraints: £,x ^ 10~^^ for 
mxc^ < 10 MeV or ^x < 10~^° s'^ for mxc^ = 100 
GeV jlj]. However, existing constraints allow ^x to be 
much larger if the photons have a wide range of energies. 

We will also consider the case of dark matter annihi- 
lation. Neglecting the evolving dumpiness of the dark 
matter field, the comoving energy deposition rate is 



= 2fmxc^ inS,f{av){l + zY 



(13) 



where the factor of two occurs because two particles an- 
nihilate in each event and (av) is the velocity-averaged 
annihilation cross section. In this case, we can define an 
effective baryon-normalized "lifetime" via 
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{av) (1 + z)^ — 

mx V "b 



(14) 



Like decay into the transparency window, annihilation 
can be thought of as a redshift-dependent decay rate. 
Obviously, annihilation injects relatively more energy at 
high redshifts. For particle decay, the free parameter is 
mxFx, which is easy to interpret. For annihilation, it is 
the less intuitive quantity {av) /mx- Reasonable values 
for neutralino annihilation are {ay) « 2 x 10^^^ cm^^ s~^ 
and mxc^ = 100 GeV (e.g., [iy|)- For illustrative pur- 
poses, we will consider a range of cross sections up to 
{av) w 2.5 X 10^^^ cm'^ s^^ (with / = 1) - somewhat 
above the upper limit provided by existing CMB data . 
(Note that our value of {av) is four times smaller than 
that used by Ref. H, because they assumed / = 0.25.) 



C. The Mean Evolution of the IGM 

We now wish to describe how these processes affect the 
IGM; our approach in this section is similar to Ref. 0. 
First consider the effects on the globally-averaged tem- 
perature and ionized fraction. Without dark matter, the 
mean IGM temperature Tk is determined by the compe- 
tition between adiabatic cooling and Compton scattering 
of CMBphotons off of the residual free electrons. We ob- 
tain [nEm 



dTk 
dt 



-2TkH{z) + (t^ „ Tk) 



2 ri2mpC^ 

3 771 fcs 



^XXh, 



(15) 



5 



where = (3TOec)/(8(TTWCMB), "cmb is the CMB en- 
ergy density at redshift z, 771 = (1 + /ne + ^i), "^2 = 
(1 + 4/hc), and /hc is the hehum fraction by number (we 
assume heUum remains neutral for simphcity) 
The mean ionized fraction evolves according to 

d^i -2- , 

— = -ax.nn + 7/2 I -^^ — \t.xx^, (16) 

where the first term describes recombinations with an 
effective coefficient a (taken from Ref. 49) and nn is 
the total mean density of ionized and neutral hydrogen 
nuclei (note that gas clumping is negligible at the red- 
shifts of interest here). The second term results from 
dark matter decay; i^ion = 13.6 eV is the hydrogen ion- 
ization threshold. Most ionizations are caused by col- 
lisions with hot secondary photo-electrons rather than 
by direct photoionization, with a total number of ion- 
izations ~ i?^/(37eV) for a photon of energy at 
Xi ~ which agrees well with our assumption that 

Xi (l~2;i)/3. Note that, by increasing x^, dark matter 
decay also makes Compton scattering more efficient in 
equation H15|) . which provides a second channel for heat 
input. This modification to the reionization history af- 
fects the CMB temperature and polarization power spec- 
tra through Thomson scattering during the dark ages. 
The WMAP experiment requires £,x ^ 10~^^ s"-"^ 001 
for energy deposition in the optically thick regime. 

A more subtle effect of the energy injection is to set 
up a global Lya background, which can couple the spin 
and kinetic temperatures through the Wouthuysen-Field 
effect. This occurs through collisional excitations in- 
duced by fast photoelectrons (see the analogous discus- 
sions in the context of X-ray heating in js^, Isil Is^). 
Each hydrogen excitation will produce a cascade of line 
photons; the Lyman-series photon produced in the cas- 
cade will then scatter through the (extremely optically 
thick) IGM. Higher Lyn photons will eventually cascade 
further, with ~ 1/3 of them being "recycled" into Lya 
photons [2^ [53. If a fraction Xa of the input energy 
eventually goes into Lya photons, the background is 

_ C TOpC^ XaixUb . . 



where is the Lya line frequency and ni, is the number 
density of baryons. Interestingly, this can induce rela- 
tively large coupling: 

1 /2 

where ^-24 = ^x/(10^^"' sec). As we shall see, this be- 
comes potentially important at lower redshifts, when col- 
lisional coupling between Ts and Tk weakens. Estimating 
the precise coupling requires a detailed examination of 
the collisional excitation processes, which we will forego 
here given our generic and crude approach. We will sim- 
ply assume that one-half of the total excitation energy 



goes into Lya photons, i.e. Xa = Xe/2 (see [s^ for a 
more comprehensive discussion). This Lya background 
also affects the IGM temperature, but the rates are com- 
pletely negligible in the models we consider [33. l36L Is^ . 



D. Fluctuations in the IGM 

Equations (|15|l . H16() . and H17() suffice to describe the 
evolution of the global 21 cm background. But a more in- 
teresting observable is how this signal fluctuates, which is 
easier to measure because the signal monopole is severely 
contaminated by foregrounds (see [3| and references 
therein) . Equation shows that these fluctuations have 
a variety of sources. How do variations in the energy de- 
position rate source 21 cm fluctuations? We can imagine 
some simple cases for perturbations to Q. The first is 
uniform deposition. This would be appropriate if, for ex- 
ample, the decays produce particles with long mean-free 
paths that then deposit energy after interacting with a 
uniform background (such as the CMB). This is proba- 
bly not realistic, because in most cases the majority of 
energy is ultimately transferred to the IGM through rel- 
atively low-energy photons or particles with short mean 
free paths. In such a scenario, the local energy injection 
rate will be proportional to the density (because the en- 
ergy deposition rate per particle is constant), even if the 
intermediate decay products have reasonably long mean 
free paths. To distinguish these scenarios, we define a 
variable 9u such that the energy deposition rate is pro- 
portional to {l + 9u5) (i.e., 0„ = for uniformly deposited 
energy). For annihilating particles, Q cx n\ oc (1 -I- (5)^. 
If the energy is deposited immediately, we would have 
0„ = 2 (assuming that (5 <C 1 and Taylor expanding); 
in the perhaps more plausible case that the intermediate 
products have long mean free paths, 0^ = 1 would be 
more appropriate. 

The fractional perturbation to the gas temperature, 
(5t, obeys (cf., ^12.. 14]) 

A5t 2 d(5 Xi ^ XiSi — Tk 
"dT ^ 3dt mhT^ 

-\^^^^Uxh[{1~ 9,^)5 + Sri (19) 

where Si = — (1 — Xi)/xi is the fractional perturbation 
to the ionized fraction. The first term describes adiabatic 
compression or expansion, the second and third describe 
fluctuations in the rate at which energy is transferred 
through Compton scattering, and the last term describes 
how dark matter decay affects the fluctuations. (Note 
that we have ignored terms of order Xi5, because the 
ionized fraction remains small in all of the models we 
consider. We have also ignored variations in Xh and Xi 
with the ionized fraction.) The corresponding equation 
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for Si is 



dt 



-axiUniSi + S + a'Sr) 



-m\^]^-^[s. + 6(1 -0^)1 (20) 



where the first term describes fluctuations in the recom- 
bination rate (which is a function of temperature; we 
have let a' = dlna/dluT/f) and the second accounts for 
dark matter decays. Note that, in the absence of extra 
energy injection, (5, « is a reasonable approximation 
at < 150 (cf. [13)7 because recombinations are slow in 
the largely neutral gas, but it must be included in the 
more general case. 

The appearance of S, St, and Si in each of these equa- 
tions implies that the different density modes source inde- 
pendent temperature and ionization modes. This intro- 
duces a non-trivial scale dependence into the calculation 
[T3L IT^ . But for a simple estimate it suffices to follow 
only the growing density mode, S (x a, which dominates 
on the large scales that may ultimately be observable 
[1^ . Writing St = 9t{z)S and Si = gi{z)S, we have 



dz 



dg^ 
dz 



9T - 2/3 



9tT^ - 9t{T-i - Ti 



K 



1 + z rjitj 

2 m'TTlpC^ f 

KxXh 



Tk{\^z)H[z) 
(1 - Qu) + 9T 



3 mkBTK^^^'^" {l + z)H{z) ' 

_9i aXinH{l+ 9i + (^'9T) 

il + z)H{z) 

mpC^ \ £,xXi (1 - Ou)+9i 
Eion J {1 + z)H{z) 



1 + z 



(21) 



(22) 



These equations make the effects of dark matter pertur- 
bations obvious. First consider gT- Adiabatic expansion 
and compression tend to drive St — > 25/3. But when 
Compton cooling is efficient (and if Si — Q), it inputs a 
constant amount of energy per particle, driving the gas 
toward isothermality {gT — > 0). Decay or annihilation 
energy that is injected in the same way (proportional to 
the local gas density) also tends to drive the gas toward 
isothermality. On the other hand, if the energy is in- 
jected uniformly, it drives f^r — > — 1 because it must be 
shared between more particles in denser regions. 

In equation H22|l. the first term holds Si constant in the 
absence of recombinations. The second term shows that 
recombinations tend to drive gi —1, because denser 
gas recombines more quickly. The dark matter decay 
term behaves similarly to its analog in equation (|21|l : 
injection at a constant rate per particle damps out per- 
turbations in Xi, while uniform injection preferentially 
decreases the ionized fraction in dense regions. 

Finally, we must consider fluctuations in the Lya back- 
ground. If the energy is injected uniformly, is constant 
and (5a = 0. If 0„ = 1, Jq, = i5. Note the contrast with St 
and Si', the coupling efficiency depends on the total Lya 
background, not the background per particle, because 
each Lya photon scatters many times. 




FIG. 1: IGM histories for long-lived dark matter. In each 



panel, the curves take ^-24 = 1, 0.1, 10 



10" 



and 0, from 



top to bottom, (a): Ionization histories, (b): Thermal histo- 
ries. Here the thin solid curve shows T-y. 



Using the baryonic power spectrum Pss{k) and inte- 
grating equations H21|) and (|22|l . we can estimate the 
fluctuation amplitude for St, Si, and Sa as a function 
of scale at any redshift. The power spectrum of the 21 
cm fluctuations is Hjl 



P2i{k,fi)^STl {(3' + ^?)Pss{k), 



where 



P' ^ P + PTgT - PHXigil (1 - Xi) + 9uf3a- 



(23) 



(24) 



For simplicity, we will average over the fi dependence 
when presenting our results. Note that, because the 21 
cm power spectrum is simply proportional to the matter 
power spectrum on large scales, we will only present re- 
sults at a single wavenumber k = 0.04 Mpc~^. The shape 
is straightforward to compute 0, , but we are con- 
cerned with the differences induced by dark matter decay 
and annihilation. Because we have used only the grow- 
ing mode, these differences are scale-independent. We 
have chosen k — 0.04 Mpc^^ because it is in the range of 
scales most accessible to observations (see iiVl|) . 



IV. RESULTS: THE 21 CM BACKGROUND 

A. Dark Matter Decay: Optically Thick Regime 

We have integrated the temperature and ionization 
equations with initial values taken from RECFAST 
at z = 350 (dark matter decay at higher redshifts does 
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FIG. 2: 21 cm signals for long-lived dark matter. Curves take 
the same parameters as in Fig.0 (a): Fluctuation amplitude 
at A: = 0.04 Mpc~^ (note that this scale is arbitrary), (h): 
Mean (sky-averaged) signal. The thin horizontal dotted line 
shows 5Tb = 0. 



not significantly affect the history). Figure ^ shows the 
resulting thermal and ionization histories for ^_24 = 
1, 0.1, 10~^, 10~^ and from top to bottom. These cor- 
respond to lifetimes tx ^ ^ x (lO^^'-lO^'') sec, assuming 
that the decaying particle makes up a large fraction of 
the dark matter. The uppermost curve would signifi- 
cantly affect the CMB and can be ruled out by WMAP 
P, 0, but the others have almost no effect on it (the 
total optical depths to electron scattering from z > 10 
are Tcs = 0.083, 0.016, 0.0060, 0.0048 and 0.0046 in these 
models). 

The most important point of Figure ^ is that, even 
if Xi remains small, Tk can still increase significantly. 
This is of course simply because the IGM is so cold at 
these redshifts: if roughly equal amounts of energy go 
toward ionization and heating, we would expect Tk ~ 
10'^(a;i/0.1) K, far above the usual temperature. For the 
larger values of , the increased efficiency of Compton 
scattering in the more highly-ionized Universe exagger- 
ates the effectiveness of heating. This aspect is invisible 
so far as the CMB is concerned, but the 21 cm back- 
ground measures it directly. 

Figure |3 shows the resulting 21 cm signals. Panel 
(a) shows the angle-averaged fluctuation amplitude at 
k — 0.04 Mpc~^ in the different decay scenarios of Fig- 
ure ^ In this (and in all figures unless otherwise speci- 
fied), we assume 9u = ^ (i-e-i the energy deposition rate 
is proportional to the local density). Panel (h) shows the 
sky-averaged 21 cm brightness temperature as a function 
of redshift. 

First consider the solid curves, which present the 



FIG. 3: As Fig. H The dotted curves take ^-24 = 0; the 
rest have ^_24 = 0.1. The solid curves show the net signal. 
The short-dashed curves set Xa = 0. The long-dashed curve 
(shown only in the top panel) assumes 9u = 0. 

standard calculation without dark matter decay (e.g., 
11, 12J). In this case, Compton scattering becomes inef- 
ficient at z ^ 200, and the gas begins to cool below T^. 
At z > 80, the density is large enough for collisions to 
drive T5 — > Tk, so the IGM becomes visible in absorp- 
tion. Below this redshift, STi, returns to zero because 
coUisional coupling becomes inefficient as the density de- 
creases. The fiuctuation amplitude behaves similarly, ex- 
cept with a peak at z ~ 50, because it is weig hted by the 
growing density mode, S cc {1 + z)~^ [lllll2||. 

The effects of dark matter decay on STf, are straightfor- 
ward and compare well to Ref . (l^ by continually heating 
the IGM, it decreases the intensity of absorption or even 
turns it into emission. In the strongest decay models, 
the extra heating is sufficient to render the IGM visible 
at z ~ 10 even without luminous sources: this is because 
the coUisional coupling rate Xc is a sensitive function of 
Tk- Nevertheless, the peak signal still occurs at much 
higher redshifts. Interestingly, the fluctuations evolve in 
a non-trivial way. At the highest redshifts, P21 is al- 
ways small. But it is weakest when dark matter decay 
is strongest, because the increased ionized fraction helps 
keep Tk ~ T^. In all of these scenarios, the rms ampli- 
tude reaches mK levels by z ~ 50. With strong heating, 
it remains large at lower redshifts even though ST^ de- 
creases, because (i) the fluctuations continue growing, 
(ii) the increase in Tk makes coUisional coupling more 
efficient |3l| . and (iii) the Lya background continues to 
contribute. 

The cases with moderately strong heating show the 
most interesting structure, because in such scenarios the 
21 cm signal changes from absorption to emission. Near 
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FIG. 4: As Fig. ^ except for energy injection in the trans- 
parency window. The bottom solid curve in each panel 
assumes no extra energy injection. The other curves take 
^_24 = 1. The dotted, dot-dashed, short-dashed, and long- 
dashed curves take noo = 1, 10~^, 10~^ and 10~^, respec- 
tively. 



FIG. 5: 21 cm signals for long-lived dark matter with energy 
injection in the transparency window. Curves take the same 
parameters as in Fig. 01 (a.)'- Fluctuation amplitude at = 
0.04 Mpc~^ (note that this scale is arbitrary), (h): Mean 
(sky-averaged) signal. 



to (but slightly before) the crossover point (at which 
Tk — T-y), P21 also goes to zero. (When JTf, = 0, over- 
dense regions have Tk > Tj and hence are visible.) This 
transition point would be a clear signature of strong heat- 
ing from some exotic process. Unfortunately, in many 
models the clearest differences occur at z < 40, when 
confusion with (rare) luminous sources may make it dif- 
ficult to separate the signal. 

Figure 13 illustrates how the different processes shape 
the curves in Figurc|21for a fiducial model with ^_24 = 0.1 
(solid curve). Here the dotted curve assumes the stan- 
dard recombination history (^-24 = 0). The short- 
dashed curve ignores the Lya photons created through 
coUisional excitation. This makes no difference at high 
redshifts, where coUisional coupling is already efficient, 
but decreases the mean signal by ~ 25% and the rms 
fluctuations by ~ 20% at low redshifts. In this regime 
the Lya background helps to maintain contact between 
Ts and Tk (even though oc [1 + z]^/^). The long- 
dashed curve assumes 9u = 0, so that the energy is de- 
posited uniformly. (Spatial fluctuations in the dark mat- 
ter decay rate obviously have no effect on STh, so we 
do not show this curve in the bottom panel.) This in- 
creases the fluctuation amplitude at higher redshifts by 
decreasing the heating rate in dense gas so that it absorbs 
more strongly. For the same reason, the "zero-point" for 
the fluctuations actually follows that of STi,. Of course, 
once the gas appears in emission, uniform energy injec- 
tion tends to damp out the fluctuations (because dense 
regions remain colder and hence less luminous). 



B. Dark Matter Decay: The Optically Thin 
Regime 

Figure 01 shows the thermal and ionization histories 
for several scenarios where energy is deposited in the 
transparency window. We take ^_24 = 1 and as- 
sume a roughly constant cross section. The dotted, 
dot-dashed, short-dashed, and long-dashed curves take 
'''100 = 1, 10~^, 10~^ and 10"'^, respectively; note that 
at any given redshift t and £,x are essentially degenerate 
in our simple model: the net energy deposition rates at 
z — 100 in these models precisely equal the corresponding 
curves in Figure ^ It is only the redshift evolution that 
changes. We note that = 0.047, 0.018, 0.0062, and 
0.0046 for these models (including only gas at z > 10). 
Thus the dotted curve substantially affects the CMB (and 
can already be ruled out), but the others have quite weak 
effects on it. 

For a fixed effective energy deposition rate, the ma- 
jor difference within the transparency window is that a 
larger fraction of the heating and ionization occurs at 
relatively high redshifts, with Xi and Tk increasing more 
slowly at lower redshifts. As a result the crossover point 
Tk > Tj occurs earlier (if it occurs at all, of course), so 
the features in the 21 cm signal (shown in Figure [S]) oc- 
cur at significantly lower frequencies. Nevertheless, they 
have the same general structure compared to the stan- 
dard calculation, with a reduced fluctuation amplitude 
at higher redshifts but stronger fluctuations at lower red- 
shift. The overall magnitude of the effect is comparable 
to the optically thick models in Figure |21 although here 
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FIG. 6: As Fig.Q except for dark matter annihilation. The 
bottom solid curve in each panel assumes no extra energy 
injection. The others take mxc^ = 100 GeV and (av) = 
1.25, 5.6, and 25 x 10"^'^ cm^ s"^ (dotted, short-dashed, and 
long-dashed curves, respectively). 



FIG. 7: 21 cm signals for annihilating dark matter. Most 
curves take the same parameters as in Fig. |^ with 9u = 1- 
The exception is the dot-dashed curve in the top panel, which 
is identical to the long-dashed curve except that it assumes 
6u = 2. (a): Fluctuation amplitude at k — 0.04 Mpc~^ (note 
that this scale is arbitrary), (b): Mean (sky-averaged) signal. 



the observable consequences tend to be stronger at higher 
redshifts and weaker at lower redshifts. 



C. Dark Matter Annihilation 

Figures and |7| show results for several models of 
dark matter annihilation, with mxc^ ~ 100 GeV and 
{av} = 1.25, 5.6, and 25x10^26 ^^^^3 g-i (from bottom to 
top) . The uppermost curve corresjDonds to the strong an- 
nihilation case considered by Ref. 8' and is near the upper 
limit set by WMAP. As expected from oc (1 -I- z)^, 
most of the energy is deposited at high redshifts. In 
Figure 03, this manifests itself as an elevated ionization 
fraction at z ^ 100, which decreases slowly toward lower 
redshift because £,^'^{z) becomes almost negligible and 
recombinations take over. The effect on Tk is much 
smaller at high redshifts because most of the energy is de- 
posited while Compton scattering is still efficient, which 
is much stronger than dark matter annihilation. On the 
other hand, at lower redshifts the gas has cooled suffi- 
ciently that even the tiny energy injection rate from an- 
nihilations suffices to increase the temperature by up to 
a factor of a few. 

Figure [7| shows that the effects on STi, and the 21 cm 
fluctuations are also relatively modest. The strong anni- 
hilation case decreases both the mean signal and the fluc- 
tuation amplitude at high redshifts by a factor of several, 
because the IGM temperature remains closer to T^. It 
modestly increases the fluctuations at lower redshifts be- 
cause of the increased Tk', at temperatures < 100 K, the 



collisional coupling efficiency is quite sensitive to temper- 
ature [sij . Unfortunately, more realistic scenarios with 
weaker (av) have considerably smaller effects that will 
be difficult to observe given the challenges posed by the 
observations; they are at best comparable to some of the 
slower decay models above. 

As with dark matter decay, perturbations in the en- 
ergy injection rate affect St and Si (and hence the 21 
cm fluctuation signal). Most of the curves in Figure [7|a 
assume 9^ = 1, or in other words that the energy in- 
jection rate is proportional to the local density. But of 
course the annihilation rate is actually proportional to 
; if the products have a short mean free path and in- 
teract with the local IGM, we would have 9u = 2. This is 
shown for (av) — 25 x 10~^^ cm'^ by the dot-dashed 
curve. Comparison to the long-dashed curve shows that 
this scenario further suppresses the high-redshift fluctua- 
tions, because it injects even more energy into the dense 
spots. Their temperatures therefore approach T^. On 
the other hand, the effects at low redshifts are negligible, 
because the IGM is approaching isothermality anyway. 

In summary, we do not expect the 21 cm background 
to be as useful for studying dark matter annihilation (as 
opposed to decay). This is because a large fraction of 
the energy is injected at high redshifts, when the 21 cm 
background vanishes. Fortunately, this is also precisely 
the regime to which the CMB is most sensitive, so to- 
gether these techniques will offer a useful window into 
dark sector processes. 



10 



V. Ha FORMATION AND IGM CHEMISTRY 

Atomic species in gas of primordial composition lack 
low-lying energy levels. Thus, before the production of 
metals, molecular hydrogen is a critical coolant vital for 
star formation to proceed. An increased H2 abundance 
could also in principle change the temperature of the 
gas and affect 21 cm observations. Computing the pro- 
duction of H2 in the early universe has a long history 
stretching back to the first calculation by Saslaw & Zipoy 
[Tfij . but more recently it has been the subject of intense 
study as computational resources and estimates of reac- 
tion cross-sections have improved. Generally, the amount 
formed in pregalactic gas, x^^ ~ 10~^ jl7, 55, 56, 57], is 
insufficient to be of physical consequence, although the 
abundance x^^ ~ 10~^-10~'^ formed in the denser and 
hotter gas of collapsed halos (where reactions can pro- 
ceed more quickly) , is sufficient to trigger star formation 
[TsL [T9I I . If dark matter decay preheats and ionizes the 
IGM, this will boost H2 formation both due to the in- 
creased temperature (since reaction rates are highly tem- 
perature dependent) and the increased abundance of free 
electrons (which are a key catalyst). It would be ex- 
ceedingly interesting if this boost was sufficient to affect 
subsequent structure formation, or left an observational 
signature. Recent estimates differ as to the magnitude 
of this boost 0, 122, 113. Here, we provide a more de- 
tailed look at H2 formation and consider some hitherto 
neglected effects. 

There are two main intermediaries by which H2 is 
produced in the gas phase: H~ and Hj. In contrast 
to previous studies, Ref. [H^ resolve all 423 rotational- 
vibrational levels of the ion and show that the 
pathway is greatly suppressed. This is because newly 
formed ions are photo-dissociated by the CMB before 
they can decay to the ground state or undergo charge 
transfer to become H2 molecules. The key formation 
pathway is therefore through H^, via the reactions: (1) 
H-l-e" ^11^+7, and (2) + H ^ H2-l-e". A net sink 
of II~ ions is mutual neutralization: (3) H" -I- H"*" ^ 211. 
The II2 production rate can be then obtained b y as sum- 
ing that the H" ion takes its equilibrium value 



y as s' 

m 



• form 



kik2XeX^in^ 



k2Xmn + k- 



k3x[W 



(25) 



Here, fci, fe, ^3 refer to the forward reaction rates of re- 
actions (1),(2),(3), fc_i is the photo-detachment rate of 
H^ due to both the CMB and nonthermal spectral dis- 
tortion photons from recombination, and n is the total 
proper density of hydrogen nuclei. We use the fits for 
^1,^2,^3 from |55j and calculate k^i as in |57| |. In par- 
ticular, we take care to compute the non-thermal spec- 
tral distortion to the CMB from HI two-photon decay 
and Lya photons from cosmological recombination (the 
two have roughly comparable contributions) by running 
RECFAST [43 and performing the appropriate integrals 
[SSl |59j. We use the parametric fit to the two-photon 
profile given by |60l | . This component generally becomes 



important at 70 < z < 120, when CMB photons above 
the H^ photo-detachment threshold energy of 0.74 eV lie 
in the steeply declining Wien tail, and before a signifi- 
cant fraction of distortion photons redshift below thresh- 
old. As we shall see, this is also the period of peak H2 
production, even in cases where the IGM is pre-heated 
and ionized. Finally, we also include coUisional H2 de- 
struction processes, with reaction rates as given in 6 J, 
generally the most significant being the charge exchange 
reaction H2 -I- H+ H^ -I- H. While normally unimpor- 
tant, they can become significant when the IGM is heated 
above ~ 3000 K. 

Before proceeding to calculate H2 abundances in our 
energy injection scenarios, we need to consider the impor- 
tance of other sources of radiation. The energy injected 
into the IGM from dark matter decay/ annihiliations in- 
troduces an additional non-thermal component to the 
radiation field, namely hydrogen Lyman series photons 
from atomic excitations (as well as two-photon decay 
products from ionized hydrogen atoms). Photodissoci- 
ation of H~ and H2 could in principle retard H2 forma- 
tion. In particular: (1) Lya (and other) photons could 
photodissociate H~ (which has a 0.74 eV threshold, and 
peaks at 1.4 eV). However, this effect is strongly sub- 
dominant to the CMB, and (at z < 120) the spectral 
distortion from cosmological recombination ai z = 1000. 
The latter is easy to see: there must be at least one non- 
thermal photon per baryon from recombination. On the 
other hand, if Xi ~ Xei then as a crude ffist approx- 
imation there are ~ <C 1 non-thermal photons per 
baryon from dark matter decay, much less than the con- 
tribution from recombination. (2) Lyman- Werner (LW) 
photons in the 11.2-13.6 eV band could photo-dissociate 
H2 molecules. Unlike H~ photo-dissociation, this possi- 
bility only exists with high-redshift energy injection: the 
non-thermal distortion from recombination only extends 
up to 10.2 eV, and rapidly redshifts to lower energies. 
On the other hand, Ly/3 and higher series photons from 
atomic excitations can photo-dissociate H2. We can es- 
timate this effect as follows. About ~ 15% of excita- 
tion energy goes into Ly/3 photons (or xp ~ fpXe ~ 5% 
of injected energy), with a much smaller fraction going 
into higher order transitions (see Each Ly/? photon 

scatters rigcat ~ 10 times before being downgraded to a 
Lya photon [s^ ]. Thus, Ly/3 photons are produced at a 
rate = mpc'^npXf3^x/{hi'i3), and destroyed at a rate 
= [{npajscj/nscatlnfj- Since the destruction timescale 
is much shorter than the Hubble time, we can consider 
the Ly/3 photons to assume their time-independent value 
and obtain n^j by setting = iip. The radiation field 
in the LW bands is then given by j'"^ — ft,c/(47r)n^, or 
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21 'mpC^XlsS.X rtscat 
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~ 10 
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= 2 X 10"'^C-24 

J/(10~^^ ergs~^ cm~^ sr" 
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(26) 



Note 
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there is no explicit redshift dependence, apart from pos- 
sible dependence of on redshift. Since fc^iss — 1-6 x 
IO'^^JI'^s"^ 62], this implies a dissociation rate per 
Hubble time of: 



fdii 



H{z) 

2.6 X 10^2 



(27) 
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Thus, at most a few percent of H2 molecules will be 
photo-dissociated by the ambient LW radiation field, and 
we will ignore this effect (of course, star formation seeded 
by H2 production could produce a much larger UV back- 
ground with much more significant consequences for H2 
chemistry) . 

Could we observe the non-thermal spectral distortion 
to the CMB produced by Lya photons from energy in- 
jection at high redshift? If Xe ~ Xi: then the comoving 
number density of Lya photons is n-y XeTib- The ob- 
served intensity today from photons produced at Zg is 
then 



4^ ni+^e) 



7.8 X 10"" 



(28) 



10- 



100 



erg s ^ cm ^ sr ^ 



at wavelengths A ~ 12([l + Ze]/100) /im. By contrast, the 
observed extragalactic background light at these wave- 



lengths is ^10 °ergs 



(see Fig. 9 of |63|). 



Observing the Lya and two-photon distortions to the 
CMB from cosmological recombination, which is brighter 



by a factor 



^([l-f Ze]/1000) ~ 10 but suffers from 



foregrounds of similar strength, is likewise still extremely 
challenging, although there have been proposals to do so 

nil. 

Since these contributions to the radiation field are 
negligible, H2 production in the case of an IGM 
heated/ionized by dark matter decays/annihilation can 
then be obtained by straightforwardly integrating equa- 
tion (|25|l with the previously calculated temperature and 
ionization histories. The results are shown in Fig. |S| We 
see that the boost in H2 production can be substantial 
in dark matter decay scenarios, with values approach- 
ing ~ 5 X lO"'*, which is comparable to the maxi- 
mal asymptotic abundance of ^ in gas cooling 
via atomic transitions from T > 10"* K (the "maximally 
heated" case of long-lived dark matter with £,-24 = 1 
tends to approach this value, but H2 is rapidly destroyed 
by charge exchange reactions as the IGM temperature 
climbs past 3000 K). The latter "freeze-out" value can be 
understood from simple timescale arguments ^61i] . Note 
that we have only computed II2 formation at the mean 
density of the IGM — it will be more efficient and rapid in 
the higher temperature and denser environments of col- 
lapsed halos. Thus, close to the maximal amount of II2 
can be formed, certainly comparable to the abundance 
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FIG. 8: H2 production in the IGM for different models of en- 
ergy injection. Top panel: long-lived dark matter, with IGM 
histories as in Fig. Q curves take C-24 = 1,0.1, 10"^ 10"^ 
and 0, from top to bottom. Note the rapid destruction of 
II2 in the top curve as the temperature climbs above 3000 
K. Bottom panel: solid lines depict energy injection in trans- 
parency window, as in Fig. |1] from top to bottom, curves 
take rioo ~ 1,0.1,10"^, and 10~^. Dashed curves are for 
dark matter annihilation, as in Fig. |S| with {av) — 25, 5.6, 
and 1.25 x lO"^'^ cm^ s~^ from top to bottom. 



XYi2 ~ lO"'* needed for molecular hydrogen to trigger 
cooling and star formation in collapsed halos (e.g., see 
E^). In particular, it exceeds the H2 abundance formed 
in low mass halos (where reactions are slower due to lower 
temperatures) catalysed by the low primordial electron 
fraction ~ 2 x 10"''. Thus, some of these dark matter 
decay scenarios may be ruled out on the grounds that 
they would naturally seed a good deal of high redshift 
star formation which would violate bounds on the ob- 
served WMAP optical depth. However, a variety of inter- 
locking feedback mechanisms are at play, so it is difficult 
to make quantitative claims without further study. 

For example, one possible caveat to the claim that 
early preheating/reionization would seed early star for- 
mation is that the entropy of the IGM will suppress gas 
accretion onto halos and instead exert a negative feed- 
back effect 65] . However, early reionization by decay- 
ing dark matter differs in one crucial respect from e arly 
reionization by stars, the scenario envisaged by Ref. l65l : 
unlike early star formation, reionization by decaying dark 
matter is accompanied by a negligible LW background. 
Ref. [6^ showed that the cores of gas accreted from a high 
entropy IGM had low enough densities that a small LW 
background would suffice to destroy any H2 formed — i.e., 
the photo-dissociation timescale was much shorter than 
the H2 cooling time. By contrast, in our present see- 
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nario the LW background is neghgible and indeed, H2 
can form and survive at the mean density of the IGM. 
Thus, as long as the gas can contract to sufficient density 
that the cooling time falls below the expansion timescale, 
the effect of the entropy floor is unimportant. Of course, 
once some star formation takes place, the LW background 
rises and H2 destruction in low density cores and the IGM 
proceeds. Therefore, detailed study is necessary to un- 
derstand if an early epoch of preheating and copious H2 
production will indeed result in extensive star formation 
in violation of WMAP optical depth bounds. 

There are a few general features worth noting about 
H2 production in these scenarios. Most of the H2 in 
all scenarios is made at z ^ 100; this is the highest 
redshift at which photo-detachment from the high- 
energy tail of CMB becomes unimportant, yet where the 
IGM gas is still sufficiently dense that reactions can pro- 
ceed rapidly. This era of peak II2 production is fairly 
independent of temperature or ionization history in the 
different scenarios. The boost in H2 production is pri- 
marily due to the increased free electron fraction; for low 
Xe, we see from equation ()25|) that xjja cx x^. (indeed, 
the peak H2 abundance can roughly be estimated from 

iCHa ~ kmizf)n{zf)XeiZf)t}i{zf), whcrC km ~ fcifc2/fc-l, 

and Zf ^ 100 is the redshift at which km peaks). The H2 
formed should not have a significant effect on the temper- 
ature of the IGM: since the II2 cooling function oc T'^ 
for T < 3000 K, the cooling time is: 

where S is the gas overdensity. This is substantially 
greater than the Hubble time and the Compton cool- 
ing time tc = 1.2 X 10^{[1 + z]/100)-'^{xe/10-^)'^ yr. 
The effects of H2 cooling are only important in dense 
virialized halos. 

Apart from the possible effect of seeding high-redshift 
star formation, there are few observable consequences of 
this large amount of early II2 formation. They could 
potentially increase fluctuations in Lya coupling (due to 
the consumption of LW photons in photo-dissociation re- 
gions), but this is likely difficult to detect. 



VI. DISCUSSION 

The 21 cm transition is, at least in principle, a win- 
dow into the dark ages of structure formation at z > 50. 
We have argued that, because of the overall simplicity of 
the (known) physics at that time - the expanding Uni- 
verse, hydrogen recombination, and linear gravitational 
growth - it presents a unique probe of both cosmology 
[llL [l3| and exotic processes such as dark matter decay 
and annihilation. The heating and ionization induced 
by the decay (or annihilation) products can significantly 
affect the IGM. These processes can modify the CMB 
power spectrum P, 0, j but this is a relatively insensi- 



tive measure because it requires a large Xi for scattering 
to be significant. 

Here we have shown that the 21 cm history is a sensi- 
tive measure of decay and annihilation during the dark 
ages (see also 0|), because it directly measures the ther- 
mal history of the IGM. In the standard calculation, adi- 
abatic cooling drives Tk to such low levels that heat- 
ing the IGM significantly requires much less energy than 
ionizing it. We have shown that dark matter with life- 
times 10^'*-10^^ sec can substantially affect the 21 cm 
background (see Figs. |21 and [SJ- These timescales are 
about three to four orders of magnitude longer than those 
probed by the CMB. The improvement is considerably 
smaller for annihilation scenarios, because those tend to 
inject most of their energy nearer the surface of last scat- 
tering, when the CMB is more sensitive and the 21 cm 
background vanishes. 

We have made predictions for both the sky-averaged 
signal - which measures the total energy deposition rate 
flSj l - and the 21 cm power spectrum. Because the tem- 
perature, ionization fraction, and Lya background all af- 
fect the 21 cm signal, the overall amplitude and redshift 
evolution (at any scale) of P21 can provide powerful con- 
straints on such exotic processes: reasonable scenarios 
produce order unity effects on the power spectrum and 
also introduce non-trivial redshift dependence if the IGM 
ever becomes hotter than the CMB. Thus any measure- 
ment of JVi during the dark ages (such as those advo- 
cated by [ll|) will be useful in this context, even if it 
only measures fluctuations on rather large scales. 

More detailed measurements can begin to constrain 
the decay and annihilation processes themselves by in- 
ferring the properties of the products that interact with 
the IGM. For example, if the decay produces either soft 
( < 3 keV) or extremely hard ( ^ 10 GeV) photons, the 
energy will be deposited into the IGM nearly instanta- 
neously 0. In this case the effects will be most obvious 
at lower redshifts, when there is more time for heating 
to take place and Compton coupling with the (spatially 
uniform) CMB is weaker. On the other hand, if the IGM 
is optically thin to the products (or if the dark matter an- 
nihilates), the consequences are more confined to higher 
redshifts (where the optical depth is larger because of the 
increased density). More subtly, the interaction processes 
between the IGM and the decay products determine the 
fluctuations in the ionizing and heating rates, which in 
turn affect P21 (see Figs. |21 and [T)). In some cases, this 
could even introduce extra spatial dependence into the 
power spectrum, although we have not examined such 
effects here. 

Rather than examining speciflc particle physics mod- 
els, we used a generic and flexible formulation for decay 
and annihilation. We refer the interested reader to the 
discussions in ^III Al and in Refs. 0, for the thermal 
and ionization histories in specific models. In general, 
however, we note that among recently popular models, 
decaying light dark matter (such as axinos and sterile 
neutrinos) would have the strongest effects. The decay 
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of heavy dark matter could also affect the thermal his- 
tory, though this depends strongly on the allowed decay 
channels. Neutralino annihilation, or the annihilation of 
light dark matter, could also provide an interesting sig- 
nal. 

We have also shown that the heat deposited in the 
IGM, as well as the excess ionization, can affect the chem- 
istry of the IGM. In particular, the increased tempera- 
ture and ionization can dramatically increase the rate of 
H2 formation. While the resulting abundance can be up 
to orders of magnitude higher than the standard value 
(and in some cases comparable to the maximal asymp- 
totic abundance xh^ ~ 10^"^ for gas phase H2 formation), 
it is not an important coolant at the low densities of the 
IGM. Thus its direct observable effects are small. How- 
ever, the increased H2 abundance will strongly simulate 
early star formation in dense halos, quite possibly vio- 
lating WMAP constraints on Tg. With better modeling, 
observations of the first stars and Te may also provide 
limits on dark matter decay and annihilation. 

We must of course acknowledge the tremendous diffi- 
culty posed by 21 cm observations at z > 50, the regime 
in which dark matter decay signatures would be cleanest. 
The principal challenge is the enormous brightness of the 
Galactic synchrotron foreground, which has a brightness 
temperature Tgky ^ 10^ K at the relevant frequencies of 
^ 30 MHz. It will make measurements of the smoothly 
varying dT^ extremely difficult; searches for fluctuations 
will probably be much easier (though still well beyond 
current capabilities). The largest transverse wavenumber 
observable by an array distributed in a circle with radius 
i?max is A:j.,max « 0.2 [50/(1 + z)](i?„iax/2km) Mpc-i. 
For a crude sensitivity estimate, we consider an array 
with uniform baseline coverage observing the spherically- 
averaged signal. Then the error on the power spectrum 
at wavenumber k would be |6(3] 



fc3(5P2i 0.1 mK / k \ ' Y Tsky 2km\ 




where /cov = ^e/ (■"'-Rmax) is the array covering factor, Ae 
is its effective area, B is the bandwidth of the observa- 
tion, tint is the total integration time, and we have binned 
the data in segments of logarithmic length ek. Here we 
have assumed that k is much larger than the wavenumber 
corresponding to the total bandwidth of the experiment 
(~ 0.015 Mpc"i for B = 10 MHz) and much smaller 
than fc_L_iiiax- Most importantly, we have assumed a large 
field of view - of order one steradian - corresponding 
to baselines of order 10 m; thus, the instrument must 
be composed of many small antennae. Our parameter 
choices in equation H30|l allow a direct comparison with 
the figures in this paper at z ~ 50 if i?max = 2 km. 
Clearly, several square kilometers of collecting area are 
required to produce any useful limits. Moving to larger 
scales can help slightly, but foreground removal will prob- 
ably compromise measurements at fc < 0.01 Mpc~^ (see 
§9.3 of Ref.l^and references therein). Fortunately, more 
compact array designs improve the sensitivity over a uni- 
form baseline distribution by factors of a few (see Fig. 6 of 
Ref. IgS). Exploration of the highly- redshifted 21 cm sky 
is just beginning, and over the next few years we should 
learn much more about what is possible. When experi- 
ments to open up the dark ages do eventually come along, 
they will provide important constraints on dark matter 
decay and annihilation - which can plausibly have order 
unity effects on the 21 cm signals. 
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